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Introduction 

For the past 20 years radiometers aboard orbiting 
satellites have provided the most comprehensive sets 
of radiation measurements of the outgoing longwave 
radiation (OLR) of the Earth on both a regional and 
a global scale. Beginning in 1975 the third generation 
of satellite experiments for measuring the radiation 
budget began. This new series of experiments was 
referred to as the Earth radiation budget (ERB) 
experiment (Smith et al. 1977 and Jacobowitz et al. 

! 1979). 

The first ERB experiment was launched into a 
nearly circular Sun-synchronous Earth orbit aboard 
the Nimbus 6 satellite on June 12, 1975. The ERB in- 
| strument package included two Earth-viewing wide- 
I field-of-view (WFOV) radiometers, one of which 
1 measured total irradiance and the other of which 

measured shortwave (SW) irradiance. The difference 
between the two (Total — SW) is the OLR. These 
WFOV radiometers viewed the entire Earth’s disk 
from the Nimbus 6 satellite, which had an altitude of 
1100 km and Equator crossing times near 12:00 a.m. 
(noon) and 12:00 p.m. (midnight). The SW channel 
had a spectral range of 0.2 to 3.8 //m, and the total 
channel measured the irradiance from 0.2 to 50-b/^m. 
The first data were taken by Nimbus 6 on July 2, 
1975. Nimbus 6 continued to collect usable Earth 
radiation data through June 1978 for a continuous 
3- year data set. 

This paper is an atlas of monthly averaged OLR 
results estimated at the top of the atmosphere (de- 
noted by TOA and defined herein to be an altitude 
of 30 km) for the 3 years of Nimbus 6 operation span- 
ning the time period from July 1975 to June 1978. 
Jacobowitz et al. (1979) analyzed and published some 
results of the first 18 months of ERB measurements 
from the Nimbus 6 WFOV data set. They used the 
inverse square approximation (geometric shape fac- 
tor) to determine the radiant exitance at the TOA. 
No formal publication of the ERB WFOV results 
of OLR from January 1977 to June 1978 has been 
released. 

The purpose of this atlas is to document all the 
WFOV OLR results from the 3 years of Nimbus 6 
operation in a form that allows analysis of the radi- 
ation field of the Earth. The results contained in 
this atlas were derived with a deconvolution (i.e., 
a resolution enhancement) technique which repre- 
sented the WFOV monthly averaged OLR as an ex- 
pansion of spherical harmonic coefficients (Smith and 
Green 1981). Tables of these coefficients along with 
monthly averaged contour maps of OLR results for 
3 years are included. Contour maps and spherical 
harmonic coefficients for the first year of the Nim- 


bus 6 ERB experiment WFOV results have previ- 
ously been documented (Bess, Green, and Smith 
1980). However, the results for the first year in this 
atlas show some improvement over earlier results be- 
cause of the use of a correction for degradation of the 
SW channel. 

The results documented in this atlas are impor- 
tant for a number of reasons. One relates to the data, 
which are both broadband and WFOV data. Because 
the measurements are broadband, the ERB radiome- 
ters offer some significant advantages over the instru- 
ments aboard NO A A operational- type polar orbiting 
satellites, which measure upwelling radiation in the 
narrow spectral regions (0.5 to 0.7 fim in the visible 
region and 10.5 to 12.5 (i m in the infrared region). 
Winston et al. (1979) published an atlas which doc- 
umented 4 years of OLR results from the scanning 
radiometer (SR) on the NOAA operational satellite. 
Janowiak et al. (1985) also published an atlas of OLR 
results derived from NOAA operational satellite data 
which covered the time period from June 1974 to 
November 1983. The primary disadvantage of mea- 
surements made by NOAA operational satellites is 
that measurements are made in narrow spectral re- 
gions which must then be empirically corrected to es- 
timate the broadband OLR. Their big advantage is 
their very high spatial resolution compared with the 
WFOV radiometer, which is limited to large scales. 

The WFOV instrument is also well suited for 
measuring large-scale features since its field of view 
is such that the instrument measures all incident 
radiation from horizon to horizon. Herein lies one 
of its advantages since the measured radiation is 
integrated over a broad variation of angles, and thus 
is less sensitive to directional models than are data 
measured with narrow-field-of-view radiometers. 

In regards to stability of instruments, scanning 
radiometers and fixed- WFOV radiometers tend to 
be very stable over time in the infrared region of 
the spectrum, but they experience degradation in 
the shortwave region. However, because of their 
mechanical simplicity, fixed- WFOV radiometers typ- 
ically have greater longevity than do scanning ra- 
diometers. For instance, the fixed- WFOV radiome- 
ter on Nimbus 6 operated uninterrupted for 3 years, 
and the Nimbus 7 fixed- WFOV radiometer has been 
operational for over 7 years. 

In addition to the data being broadband and 
WFOV, the method of representing the data using 
spherical harmonics is important. The spherical har- 
monic coefficient data set for each month represents 
a condensation of the OLR field. Some of these co- 
efficients have a physical interpretation attached to 
them. These coefficients may be analyzed individ- 
ually or in combination to study different aspects 


I 



of the radiation field. The results from Nimbus 6 
also represent the beginning of a long-term time se- 
ries of OLR for the ERB WFOV-type instruments. 
This 3-year data base takes us into the time period 
of the second ERB WFOV radiometer, which was 
flown on the Nimbus 7 satellite in November 1978. 
The Nimbus 7 data are the subject of another atlas 
of OLR. Nimbus 7 has been operating continuously 
for 7 years, and when its data are combined with the 
data base from Nimbus 6, we have a time series of 
OLR results covering 10 years. This time series data 
will be very valuable for doing monthly, annual, and 
interannual studies of OLR. 

The authors are grateful to H. Jacobowitz of 
NO A A for providing the Nimbus 6 data tapes which 
made this research possible. Also, many thanks to 
T. P. Charlock, J. J. Buglia, and W. F. Staylor for 
their many suggestions and to M. A. Woerner for her 
programming help along the way. 

Data Processing and Analysis 

The 3 years of ERB WFOV irradiance data from 
the Nimbus 6 satellite covering the time period from 
July 1975 to June 1978 were supplied by NOAA. 
Only the OLR data are documented in this atlas. 

The nominal duty cycle of the ERB radiometer 
was 2 days on and 2 days off. When the radiome- 
ter was turned on, measurements were taken at 4- 
sec intervals along the orbital track. The data tapes 
supplied by NOAA contain data reduced by averag- 
ing four consecutive measurements, resulting in one 
averaged value every 16 sec. 

Because the Nimbus 6 ERB radiometer was con- 
strained to operate on a duty cycle of 2 days on fol- 
lowed by 2 days off, continuous daily measurements 
were not possible. This sampling strategy tends to 
change the true radiant exitance. However, the ab- 
sence of a continuous daily data set did not seriously 
constrain monthly averaged WFOV data since the 
data were smoothed in the averaging process over 
1 month and over the large spatial area of the WFOV 
radiometer. Green and Smith (1978) looked at the 
temporal variation over six duty cycles of 2 days each 
for 1 month of Nimbus 6 data. Their results showed 
very little change from one duty cycle to another on 
a global and zonal scale. Small changes occurred for 
some regions. 

The other inherent sampling bias characteristic 
of all Sun-synchronous polar orbiting satellites is 
that they measure OLR at only two local times. 
For Nimbus 6, 12:00 a.m. (noon) and 12:00 p.m. 
(midnight) were the equator crossing times. Because 
of this sampling bias, diurnal variations cannot be 
studied. 


These tapes were processed by taking daily mea- 
surements and averaging over 1 month and over 5° 
increments in latitude that formed an igloo-type grid 
system of near-equal-area regions. The igloo grid 
system is symmetrical about the Equator and has 
3 grids at the polar regions from 85° to 90° latitude 
and 72 grids at the equatorial region from 0° to 5° 
latitude. Figure 1 is a sketch of the igloo grid system 
for the Northern Hemisphere showing the number of 
grids in each 5° latitude region. The total number of 
near-equal-area grids is 827. The igloo grid system 
is symmetrical about the Equator, giving a total of 
1654 grids over the Earth. 



Figure 1. Igloo grid system for near-equal-area grids. Num- 
ber of grids in each 5° latitude region for one hemisphere 
is shown. 


Some corrections and editing had to be made 
to the ERB WFOV measurements before they were 
suitable for data analysis. Sun contamination is a 
problem which occurs near sunrise and sunset with 
the WFOV radiometer because the field of view is 
larger than the Earth’s disk. Measurements were 
eliminated when the solar zenith angle at the nadir 
(or subsatellite) point was between 99° and 123°. 
This was the range of angles deleted by investigators 
when analyzing Nimbus 7 ERB WFOV data (Kyle 
et al. 1984). When the first year of Nimbus 6 
WFOV data was analyzed (Bess, Green, and Smith 
1981), solar zenith deletion angles were from 111.5° 
to 123.5°. However, with the range of solar zenith 
angles specified by the Nimbus 7 investigators, it was 
found that another possible source of error in OLR 
measurements, caused by a thermal transient in SW 
measurements because of irradiance of the sensor at 
spacecraft sunset, was eliminated. Other editing had 
to be performed occasionally because of anomalies 
in the data that have been well documented (Bess, 
Green, and Smith 1981). 

In addition to the expected editing which had 
to be made to the WFOV measurements, the Nim- 
bus 6 ERB WFOV radiometer had a systematic error 
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which caused the WFOV measurements to be low by 
about 11 percent (Smith et al. 1977). This was the 
case for both the total and the reflected solar radi- 
ation channel. As such, all measurements from the 
tapes at satellite altitude were increased by 11 per- 
cent as a calibration correction. 

A correction was also applied to the measure- 
ments to account for degradation over time of the SW 
radiometer. The degradation appears to be strictly a 
SW radiometer problem, since no measurable degra- 
dation occurred during the nighttime portion of the 
orbit of the total channel radiometer. It seems to be 
characteristic of S W radiometers to experience degra- 
dation over time, not only on ERB satellites, but also 
on NOAA operational satellites. The degradation, 
which was confined to the ascending portion of the 
orbit, caused the OLR to increase in a linear fashion 
over the 36 months for which Nimbus 6 WFOV mea- 
surements are available. The correction for degrada- 
tion was different for different latitude zones, being 
generally larger in the tropics and in the middle lati- 
tudes. The net result of degradation was to cause the 
OLR measurements to increase over time at a rate of 
about 3 to 4 W/m 2 per year. 

To account for the degradation, a ratio-to- 
centered moving average method was used for 5° lati- 
tude regions to remove seasonal effects from the data 
prior to performing any trend analysis (Smith and 
Williams 1971). This was necessary so that seasonal 
variation could be removed before removing the trend 
from the data. Basically, the ratio-to-centered mov- 
ing average method determines a monthly seasonal 
index by dividing the original data for a given month 
by the 12-month centered moving average and tak- 
ing the average of that month over time. This was 
done for each of the 12 months of the year. Dividing 
each value of the original data set by the appropriate 
seasonal index gives the percentage by which a given 
monthly average was above or below the average for 
the year. After the set of data has been adjusted for 
seasonal variation, the resulting data give a better 
representation of the trend and the random effects. 

After the data were adjusted for seasonal varia- 
tions, linear least-squares procedures were applied to 
the data to determine the slope of the trend line for 
the time scries in each latitude region. The slope 
gives the rate of degradation in watts per square me- 
ter per month in the OLR and is used to determine 
the correction in OLR for a given month and latitude 
region. Figure 2 shows these results as a plot of the 
rate of degradation per month against latitude. The 
rate is greatest from 60°S to 30°N latitude and peaks 
between the Equator and 30° N latitude. The rate of 
degradation drops off at the higher latitudes. 



Figure 2. Degradation of shortwave spectral region per 

month. 

After the measurement data were corrected for er- 
rors and were time- and space-averaged to obtain re- 
gional monthly averages, a deconvolution (resolution 
enhancement) technique was applied to represent the 
radiant exitance at the TOA by a truncated series of 
spherical harmonics. The deconvolution technique 
takes advantage of the fact that spherical harmonics 
are the eigenfunctions of the measurements opera- 
tor and reduces the radiant exitance field from satel- 
lite altitude to TOA by dividing by the appropriate 
eigenvalues. All the results in this atlas are based on 
this deconvolution technique (Smith and Green 1981 
and Bess, Green, and Smith 1981). 

The governing equation from which the monthly 
spherical harmonic coefficients were produced is 

N n 

M(0,<M) = H N nPn(cos9) [<C(«)coe(m*) 
n=0 m=0 

+ S™ (t) sin(m0) ] 

where M(6 , 0, t) is the measurement at the TOA, 9 
is the colatitude, 0 is the longitude, t is the time, and 
P™ is the associated Legendre polynomial of degree 
n and order m. The terms C™(t) and S™(t) are 
the even and odd real spherical harmonic coefficients, 
and the normalizing factor is 

N™ = [(2 n + 1 )(n - m)! (2 - Stf) f(n + m)!] 1 / 2 
where 6™ is the Kronecker delta function. 

Discussion of Results 

Included in this atlas are spherical harmonic 
coefficients and associated global contour maps of 
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36 months of Nimbus 6 ERB WFOV outgoing long- 
wave radiation results. It is not the intent in this at- 
las to do an in-depth analysis of the data but rather 
to compile and document the spherical harmonic co- 
efficients and the associated contour maps which can 
then be used to do many kinds of valuable analysis. 

Each table in this atlas contains a set of spher- 
ical harmonic coefficients for 1 month of mean val- 
ues. Results are for a spherical harmonic expansion 
truncated to the 12th degree. For such a 12th-degree 
expansion, 169 coefficients are required to specify the 
radiation field. The coefficients above the stair-step 
line are the 78 sine terms. With the exception of the 
first column, the coefficients below the stair-step line 
are the 78 cosine terms. The first column contains 
the 13 zonal terms. The sine and cosine terms repre- 
sent the nonaxisymmetric terms and give a measure 
of longitudinal variation. The format of the tables 
makes it very easy to pick off any coefficient. The 
superscript m is the longitudinal wave number or 
order and n represents the degree of the spherical 
harmonic. Thus, in the first column, which repre- 
sents the zonal terms, m is 0 and n ranges from 0 to 
12. Physical interpretations can be associated with 
some of the zonal terms. Thus, Cq is the global aver- 
age, C? is a measure of hemispherical or pole-to-pole 
difference, and C® is a measure of Equator-to-pole 
gradient. It has been shown that over 80 percent of 
the degree variance is in the zonal terms (Smith and 
Bess 1983). This variance is because at large scales, 
Earth-emitted radiation is strongly dependent on 
latitude. 

The monthly averaged spherical harmonic data 
sets can be used in a variety of ways to study the OLR 
on regional, zonal, and global scales in the spatial 
domain and on monthly, annual, and interannual 
scales in the time domain. One application is to 
model the global radiation field. The advantage of 
such models is that they can represent large data sets 
with relatively few parameters. Another advantage 
of a spherical harmonic representation is that it 
provides a mathematical structure that permits one 
to study separately the latitudinal variations using 
the zonal coefficients and the longitudinal variations 
and wave properties using tesseral coefficients. In 
short, spherical harmonic representation allows the 
radiation field to be broken into its component parts, 
which can then be studied separately or in various 
combinations. For example, the coefficients are well 
suited for time series analysis, spatial spectra studies, 
and parameterization studies (Smith and Bess 1983 
and Short et al. 1984). 

In this atlas each spherical harmonic coefficient 
set has a companion monthly averaged global contour 


map of OLR. The OLR less than 240 W /m 2 is shown 
as dotted contour lines. The contour interval is 10 
W/m 2 , and highs and lows are shown. These contour 
maps give a “quick look” of how the OLR is varying 
over monthly, annual, and interannual time scales. 
With their associated sets of harmonic coefficients, 
they form a valuable data set for studying many 
different aspects of our changing climate. 


NASA Langley Research Center 

Hampton, VA 23665-5225 

July 7, 1987 
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